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Abstract 

We study several examples of kinetically constrained lattice models using dynamically accessible 
volume as an order parameter. Thereby we identify two distinct regimes exhibiting dynamical 
slowing, with a sharp threshold between them. These regimes are identified both by a new response 
function in dynamically available volume, as well as directly in the dynamics. Results for the self- 
diffusion constant in terms of the connected hole density are presented, and some evidence is given 
for scaling in the limit of dynamical arrest. 
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In nature, molecules, particles or other elements of a nearly arrested system may some- 
times stop moving without any accompanying sharp change in the thermodynamic quanti- 
ties. Such systems have not crystallized, the free energy has not been minimized, and there 
is no obvious order parameter for the 'transition'. This phenomenon of precipitous, and 
apparently collective, loss of motion has been named 'dynamical arrest' jl} Q, Q, Q m or 
jamming p, [l,!^- It is involved in the 'glass transition' P and even, apparently, gellation 

, very significant advances have been made in 
developing our understanding 3, Q|- in the scientific community there is an emerging 
opinion that the many modes by which complex condensed states are formed are all aspects 
of the process of dynamical arrest. 

We have sought a fundamental theoretical approach to dynamical arrest based on 

a physical order parameter to describe the vicinity of dynamical arrest, keeping in mind the 
potential to make direct contact with experiment. In particular, the advent of new experi- 
ments an d new experimental methods being developed in colloidal and 
soft matter science, opens up the possibility to directly connect observation of particle con- 
figurations and their fluctuations to suitably framed theoretical concepts. We have identified 
dynamically available volume (DAV) Q], the ensemble of physical space available to particle 
dynamics, as a possible order parameter. Its use has been illustrated in a recent study of 
certain lattice models that exhibit dynamical arrest (3, 0, Q] and has been considered in 
models of granular matter 

In this paper, we highlight a new insight that for the simple model systems we explore, 
and possibly more generally in nature, the geometry of available empty space of a nearly 
arrested system has a profound impact on the nature of the relaxation processes and the 
accompanying laws governing the onset of arrest. We propose that the onset of arrest involves 
two qualitatively different regimes, both of which appear to exhibit a degree of universality 
and scaling not previously appreciated. Our observations are valid for a variety of lattice 
types, and in both two and three spatial dimensions. 

The models we study include kinetically constrained lattice models introduced by Kob 
and Andersen j3] in which particles may move to an adjacent empty site (a hole) if sur- 
rounded by c or fewer neighbors, and if the movement is reversible according to these same 
rules. Vacancies are empty sites into which no adjacent particle can move. A dynami- 
cal arrest transition occurs as a function of increasing particle density, system size, spatial 
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dimensionality d, and constraints c. 

At moderate densities a particle that moves into a hole typically generates a new hole 
that may be used for a subsequent motion. By such sequential motions, all particles in 
the system may be moved beginning with such a ('connected') hole. At high density many 
('disconnected') holes are caged by a continuous boundary of particles that cannot be broken 
by particle rearrangements inside the boundary (eg. in the square-lattice KA model the 
boundary has at least a double row of particles of rectangular shape). Each particle on 
this boundary is blocked because it cages, and is caged by, its neighbors. In our previous 
work we deemed disconnected holes irrelevant for the long length and time scale transport 
coefficients, and the connected holes (of density u r ) become the natural order parameter 

□ n 

|2jj. Subsequently |3j it has become apparent that, asymptotically close to final arrest (in 
the lattice models, for p = 1) this is equivalent, via the relation v c = to relating the 

transport coefficients to the bootstrap length £ jsjIQ- 

Connected and disconnected hole densities can be estimated by numerical simulation, 
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or by theoretical approaches j^- in some two dimensional models they can be calculated 
exactly [19|- The total hole density v t is a sum of the connected v c and disconnected Vd 
hole densities, and v c is plotted for c = d = 2 KA in Fig. Q (c = d = 3 KA in the 
inset). Although the results shown here for the connected hole densities, or equivalently, 
the bootstrap length £ d = \jv c correspond to bootstrap simulations for system sizes beyond 
any previously computed, it is important to note that they remain far too small P, Q] to 
be described by simple (exact) asymptotic approximations 2^1 . 

The connected hole density is the order parameter of the dynamical system, and deter- 
mines the transport coefficients. The response of the dynamical system to creation of a 
hole may be calculated from the derivative (essentially a response function for hole creation) 
Xc — di/ c /dv t , plotted as a function of density for the c = d = 2 KA model in Fig. |21 The 
result is striking- we find a threshold ('transition') density below which new holes automat- 
ically become connected (and therefore contribute to diffusion) and above which nearly all 
holes are trapped within extended cages. 

To the low density ('unstable') side of the transition a new hole is typically close to an 
existing connected hole, and therefore likely also to be connected. More holes lead to further 
mobility, and the system is termed unstable. However, this mechanism also implies that 
dynamical connectedness, formerly defined only in relation to phase space (the possibility 
to move all particles from a given initial condition) is strongly associated with geometrical 
aspects of a network of connected holes in this regime. 

To the high density (' met ast able') side of the transition, creation of a new hole typically 
leads to a rattler, or disconnected hole due to caging. Existing connected holes are now 
dispersed sufficiently far apart that they provide no assistance to a new hole. 

The small peak in Xc (at u*) is considered the onset of this transformation of phase space, 
and the peak of the second derivative, t] c = d 2 v c /dv^, an effective threshold at which the 
system may be deemed to have made a 'transition' in dynamics (see inset to Fig. |2j). In that 
inset we show that a similar reorganization of phase space is observable in models with dif- 
ferent sets of kinetic constraints, symmetry of lattice, and spatial dimensionality. It should 
be emphasized that it has so far been difficult to quantitatively relate this restructuring of 
phase space, and the change in dynamical processes, to simple static concepts of connected- 
ness of vacancy and hole networks, though attempts continue [24J . On the contrary, for each 
of the different models described in this inset there are quite different spatially connected 
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FIG. 2: Xc = du c /du t for d = 2 KA model. Also sho^0n are —dv c /dp (o) and —dv t /dp (A). The 
inset shows ry c = d 2 u c /duf for cubic KA c = 3 (o), fee KA c = 6 (A) and KA c = d = 2 (•). 



networks, whilst the dynamical phenomenon is quite general. Note carefully that the con- 
nected hole density contains information on the potential temporal evolution of the system. 
We now consider the important step in understanding the geometrical aspects involves the 
concept of response functions. These are defined one level removed from the specific rules 
of the model itself, that capture the spatio-temporal connected-ness of phase space. This 
approach is expected to be more appropriate also in attempts to generalize to the continuum. 

We conclude there are likely two quite distinct dynamically slowed regimes on approach to 
arrest, possessing qualitatively different spatio-temporal or dynamical relaxation processes, 
different laws for the diffusion constant, and non-linear response functions. There is ample 
anecdotal experimental information for differing dynamical regimes, and possibly differing 
types of law, but hitherto there has been no clear formulation of the underlying ideas. We 
are in a position to check this explicitly in the dynamics of these models. Note however, 
that the strength of the calculations above are that they are (or may be made) essentially 
exact, whereas the dynamics calculations are lengthy with more limited reliability. We shall 
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now establish the existence of two different regimes for these models by exploring different 
aspects of the dynamics. 

The mean squared distances travelled by particles (r(t) 2 ) are computed as usual and 
the diffusion constants D calculated from the long-time diffusivity (r(t) 2 ) = 2dDt (see Fig. 
OJ). Having computed the relation between the connected hole density and particle density 
we may eliminate the latter and plot the diffusion constant in terms of the order parameter, 
u c . It is known that the connected hole density fails to reach a simple asymptotic form in 
any density range accessible to computers 0, Q, 25, 26, 27, 28^ Indeed, in the regimes we 
discuss, there is a difference of several orders of magnitude from the asymptotic form that is 
usually assumed. Therefore, in writing the diffusion constant in terms of the connected hole 
density (or dynamical correlation or bootstrap length) it is important to use results that are 
appropriate to the density regime under consideration. 

For moderate density we recover the known [2j quadratic law D = 71/?, (7 is a model- 
dependent effective rate constant) for different models, rules sets and spatial dimensions. In 
that regime v c and p (and therefore D and p) are connected by a power law. However, at 
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higher densities, for those models that have been explored beyond the 'transition', there is 
a sharp change of slope and a new law for the diffusion constant, D = 71/! (z non-integer). 
In that regime the hole density (and thereby diffusion constant) is connected to particle 
density via an exponential-like decay, though the detailed form is subtle In each case we 
have marked that value of particle density at which the unstable-to-metastable transition 
takes place based on available volumes. Examples are given for two and three-dimensional 
KA models, but similar phenomena are present in several other lattice types, and kinetic 
constraint sets. For the three dimensional cubic lattice the cross-over is on the edge of the 
presently accessible length and time scales. 

One should not be complacent about quantitative treatment of this high density limit as, 
for example, expressed in the (effective) exponents z, which may contain a slowly varying 
density dependence. The fact that the diffusion constants are derived from some of the 
longest simulations and large system sizes is no reassurance; the dangers of interpretation 
associated with the (related) subtle asymptotic phenomena in the bootstrap percolation 
problem where system sizes far beyond current computation fail to reach asymptotic laws is 
sufficient warning 0, 0,13, 3- Here we will seek to clarify the distinct physical relaxation 
processes. The true asymptotics may have to wait for developments in dynamics that mirror 
those that have taken place recently in the bootstrap problem. 

Amongst the most striking way to illustrate the consequences in dynamics of the transition 
described in Fig. is to represent the spatio-temporal processes themselves visually. Thus 
in Fig. H]we illustrate qualitatively different dynamics of unstable and metastable regimes 
from several models using two representative densities on either side of the geometrical 
transition. 

The examples correspond to c = d — 2KA and KA Modified and c = 6, fee KA models 
(the c = 2 KA Modified model is identical to the c = 2 KA model with the added restriction 
that when we consider a move, any of the two vacant neighbours of the particle must be 
second neighbours to each-other). In each case the particles which have moved after some 
time are shown in their initial positions. We have explored many examples, and found the 
phenomena we describe to be quite general. In the unstable regime motion spreads rapidly in 
a concerted manner, with ample pathways for configurational relaxation arising from implied 
networks of connected holes. At higher particle densities, the available volume disconnects 
('de-percolates'), leading to the metastable regime. There the pathways are quite different, 
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(a) KA c = 2, p = 0.82, 0.94 (6) Modifed c = 2, p = 0.80, (c) KA FCC c = 6, p = 0.75, 

0.87 0.85 
FIG. 4: (Color online) Sample configurations of simulations of different models in the unstable 
regimes (upper panels) and metastable regimes (lower panels). We show only the particles that 
have moved after some time. The patterns in the upper panels develop almost immediately, while 
those on the bottom are essentially unchanged after many millions of MCS. These give a pictorial 
representation of how dynamically accessible volume is delocalized in the system, and thereby 
characterize the nature of the relaxation processes in the two regimes. In the unstable regime 
spinodal-like waves spread throughout the system, whilst in the metastable system movement 
initiates in localized droplets, and spreads slowly from those (rare) seeds. 

the motion is much slower and mediated by 'droplets' within which connected holes mobilize 
the particles with the assistance of a network of vacancies. These pictures, illustrating 
the mechanisms by which dynamically available volume is delocalized and configurational 
relaxation occurs are quite consistent with the understanding of the two regimes arising 
from the susceptibility (x c ) alone, affirming the value and link between dynamics, and static 
averages in available volumes, based on bootstrap processes. 

We end by noting that the metastable regime appears to exhibit a characteristic time, 
and a degree of scaling, on approach to true dynamical arrest. Thus, in Fig. El we show 
the scaled curves, and in the inset the equivalent crossover time (r) at which the system 
becomes diffusive. The scaling behaviour, and the requirement for appropriate short and 
long-time behaviour leads us to write (r 2 ) = t 2/3 /(t/r) with f(x) ~ const (small x) and 
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FIG. 5: Scaling of the mean square displacement (r 2 ) = t 2/3 /(t/r) for 2D KA in the metastable 
regime. The densities are p = 0.850(0), p = 0.860(+), p = 0.870(d), p = 0.900(x), p = 0.930(a) 
at L = 512 and p = 0.950(o) at L = 1024. In the inset we show the diffusion constant against 
crossover time. 



f(x) ~ x a (large x). Here 13 is an exponent reflecting the sub-ergodic nature of the system 



up to the cross-over time 3CJ. The two asymptotic behaviours match at the crossover time, 
leading to R 2 = Dr, R 2 = t 213 , and therefore D ~ r 2/3_1 . For large density the connected 
hole density is written in terms of the bootstrap length, and therefore we finally obtain the 
law, D ~ with /z = —{2(3 — 1)/ (3d. From the cross-over data we estimate the exponent 
/i = 1.7, in good agreement with direct measurements, and in accord with recent studies of 
a similar model, the 2 vacancy assisted triangular lattice gas [^l]. The uncertainties of using 
cross-over data are considerably larger than the direct measurement of the exponent from 
the diffusion constant data above in Fig. El but the routes are entirely independent, and 
therefore valuable in the context of the present uncertainties about finite size dependence. 
We conclude that it is possible that the cascade of cages and partial cages determining 
the subergodic behaviour provides a set of scalable 'traps', and thereby a collapse of the 
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dynamics data as outlined above. 

In summary, we propose the existence of two regimes of near-arrested dynamics with 
a transition between them driven by the underlying geometry of available volume. The 
geometrical properties of this empty space are not simple, but are characterized by a new 
response function. The two regimes possess different relaxation processes, leading to different 
dependencies of the diffusion constant on density. That change in the diffusion constants 
could for some cases, be so dramatic to be mistaken for the dynamical arrest itself. One 
would expect to be able to observe the geometrical transition in direct imaging experiments, 
and its consequences via new experimental methods that can probe long-time and long 
length-scale dynamics in colloidal and nano-particle science. 
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